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Informal Technical Information 


1. Actuator Testing 

The test results for our current actuators in the low frequency range have been 
reported in prior reports. We have recently conducted several new tests to our curr 
Terfenol actuators to study their frequency response, especially in the higher frequency 
range Using a Hewlett Packard combined impedance/spectrum analyzer, we measured 
acceleration response of the current actuators. We also monitored the impedance of the 
actuator during the test. Since the Terfenol material is a nonlinear material and it may 
exhibit different behavior at different levels of control current, we tested the trans 
StaSm at 15 volt, 30 volt, and 60 volt input si 6 nal levels. At the lower 
excitation input level, the transfer functions looked very similar. Some of results sho_ 
in Figure 1.1 There is no significant structural mode at the frequency range below 3000 

Hz. 

The prototype of our new design of actuator has also been tested. Figure 1.2 shows 
the displacement response of the prototype actuator at different input current levels, i is ^well 
f thei? cower spectrums. As shown in the figures, the response is fairly linear at lower 
current drive levels. The responses degrade with larger amplitude inputs. The test input 
signal is a 25 Hz sinusoidal. The displacement measurement device used was an LVDT. 
Since the LVDT device is essentially a DC measurement device, it own reson ^™^ e ^ 
about 10 Hz, and may effect our results. We are investigating this issue now We plan to 
use a capacitance probe to obtain more accurate measurements, and the National Institute of 
Standards and Technology has loaned us a capacitance probe for this purpose. The results 
will be reported in the next progress period. 

2. Mechanical Amplifier Design 

We have finalized our mechanical amplifier design, as shown in Figure 2. L The 
simulation data for the amplification ratio of an ideal design model is shown m Figure 2.2 
where the ratio Rra is approximately 4.5. The actual amphfication ration for our amplifier is 
expected to be 2.5 to 3.0. We have contracted with ComSat to fabneate the hardware for 

the amplifier. 


Some of design issues of the flexure based mechanical amplifier are discussed as 
follows: 


Advantages of flexures : 


( 1 ) 

( 2 ) 


There is no backlash which is usually a problem of using conventional joints for 
precision devices. 


They are wear free because there are no 
mechanisms will be fretting and corrosion at 
stresses are applied. 


sliding pairs. The only likely wear 
fabricated interfaces where changing 
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(3) They can be manufactured from a single piece of material to provide a monolithic 
structure; 

(4) Displacements are smooth and continuous at all levels. 

Disadvantages: 

( 1 ) 


( 2 ) 

( 3 ) 


There will be hysteresis due to dislocation movement in most materials. Its 
Xiwd^epends upon the stress level, temperature and the gram structure and 

atomic bonding in the material. 

They can only facilitate small displacement for a given size and stiffness. 

They cannot tolerate large loads. Unless precautions are ta^n acdden^ overloads 
can lead to fatigue, work hardening and eventually to catastrophic failure. 


Stiffness of Flexures 

A flexure joint can be modelled as a simple cantilever beam: 





The stiffness of the flexure under transverse load W can be expressed as: 


X = 


1 2 E I _ Ebd 3 


L 3 


L 3 


whprr- E is the Young's Modulus of the material; ^ 

I is the second moment of the cross-sectional area. For rectangular cross-sectio , 

I = lui! 

12 

Hinge Bending Angle Estimation 

Based on the beam theory, the hinge bending angle can be approximated as 


e = 


WL 3 
2 E I 
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where 9 is the bending angle at the tip of the hinge in radians, W is transversal load. 
Rearranging the above equation, we have : 

t y _ 2 E 1 9 _ E b d 3 ,9. 

L 3 6 L 3 

Numerical example: 

In the joint we designed 

E = 136500 (=210000*0.65) N/mm 2 ( for Beryllium copper ) 
b = 0.5 in = 12.7 mm 
d = 0.01 in = 0.254 mm 
L = 0.1 in = 2.54 mm 

0 = 1 degree = 0.01745 rad 

W = 136500 * 12.7 * 0.254 3 * 0.01745 1 ( 6 * 2.54 3 ) = 5.04 N 


Maximum Deflection Range 


The deflection range of a flexure stage will be governed either by its stiffness, 
through a maximum force that can be applied to it, or by its elastic properties. Using simple 
beam theory, this limit will occur when a maximum allowable tensile stress is reached at the 
section of maximum bending moment The maximum permissible displacement is 

_ O’max L 2 
Mmax — 


3 E d 


Referring to the above equations, increasing the thickness of the 
blade d will rapidly increase overall stiffness with little weight penalty. However the 
deflection that can be obtained without exceeding a specified maximum stress decreases 
linearly with the thickness. The width of the blade linearly effects both stiffness and weight 
without influencing maximum deflection. 


Euler Buckling Load 

The buckling should also be considered in the flexure design. 



L 


The critical load or buckling load is derived as: 

Pe = tc 2EI = E d b 3 

L 2 12 L 2 

Numerical example: 
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In the joint we designed N 

E = 136500 (=210000*0.65) N/mm 2 ( for Beryllium copper ) 
b = 0.5 in = 12.7 mm 
d = 0.01 in =0.254 mm 
L = 0.1 in = 2.54 mm 

P E = ( K 2 /12 ) ( 136500 * 12.7 * 0.2543 / 2.54 2 ) = 3621.5 N 


Material Selection 

A set of functional requirement can be established to guide the selection of 
materials. The flexure joint must have a good repeatability, implying a material with a good 
dimensional stability and low hysteresis loss. The material must be able to sustain 
considerable stress to allow a large degree of flexure. High resonance will be desirable : m 
order to reduce the control-structure interaction, so low mass and high mechanism stiffness 

are wanted. 

Beryllium copper has very good yield strength therefore allows greater deflection of 
the flexure? We are fabricating the prototype mechanic! amplifier out of benllium copper, 
however in a final design titanium would be a better choice. 


3. Neural Network Control System Development and Experimental Evaluation 

Vibration control of structures can be achieved by passive or active means. A typical 
passive damping technique, which is effective in damping out higher frequency vibration, exploits 
vfscoelas^maferials. Better performance at lower frequencies with less weight can be achieved 
us in<r active controls. However, design of an active vibration control system is a challenging task 
because of the lar^e number of structural modes involved in the bandwidth of interest, and the 
£ Sioum of Uncertainty involved. In our application, the problem is exacerbated by. die fact 
that the dvnamic behavior of these systems will hardly be known precisely during design, and 

identification of the system model will provide only limited ^SseTconc^m as 

such as dynamic perturbations, transient thermal states, parameter drift, etc. This rais 

to robustness and stability. 

There have been a number of active vibration isolation approaches proposed m recent 
vears m conffol del "n of most of these approaches relies on model information to predret die 
Sure's dynamic response. In reality, a perfect model of a lightly damped space struemre ha, 
an infinite number of lightly damped modes. Consequently, higher order modes are truncated fro 
the mode! in order to perform effective control design. The model-based control design such as 
So^pUmSl control, pole-zero cancellation, etc., will work well as long as the modes t included 1 in 
the model are accurately compensated, and the unmodelled modes are sufficiently damped. In 
case where there are closely spaced higher order unmodelled modes, the performance of 
control design will severely deteriorate, even resulting in instability. 


3.1 Neural network structure and backpropagation algorithm 
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A muitii-iver feedforward neural network model consists of an input layer, one or more 

Xk = [ xo, X lk , X2k, Xnk 1 is input vector, and the subscript k denote the adaptation index. 

W k = [w 0 k, wik, w 2 k w nkl is the weight vector which is to be adjusted during the 

adaptation to minimize the error. 



The output of the neuron is calculated by: 

y k = sgm(Sk) = sgm( W k X k ) (3,1) 

where sgm represents a sigmoid mapping function of the linear output signal s. 

The leamine nrocess of the network begins by presenting an input vector X to the network 

square of the sigmoid error defined as: 

e k = * yk = d k " sgm(Sk) (3.2) 

The instantaneous gradient estimate obtained at the kth adaptation cycle is given by. 


3 wi 


3w k 

3 s k 


3 w u 


= - 2 e k sgm'(s k ) — - = - 2 e k sgm'(s k ) X k 

3w k 

therefore for a general sigmoid function, the adaptation law is 


(3.3) 


W k+ i = W k + *i(- VwJ 

= \V k + 2|i e k sgm'(s k )X k 


(3.4) 
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where [L is a positive-value learning rate constant. 

For the purpose of error correction the output error at the output layer tteuron should be 
back-propagated into hiddett layers based on chain rules. From a single neuron po.nt of vtew, the 
error back-propagation algorithm is derived as. 


v x = 3 -i^= 2e k ^=2e k[ -i^] 

‘ b Xk 3Xa ax k 

= - 2 e k sgm’(Sk) — L = - 2 e w sgm'(s k ) W k 
5X k 

In the case where the sigmoid function is the hyperbolic tangent function tanh(.) : 

y k = sgm(s k ) = tanh(s k ) = ( 1+ *_ 2 — > 
we then have 

sgm ' (Sk) = g («anh<Sfc» _ j . (tanh(Sk)) 2 = i . (sgmfsit)) 2 = 1 - y| 
3s k 


(3.5) 


(3.6) 


(3.7) 


The adaptation law becomes: 


W k +i = W k + 2 |i e k (1 - y k ) X k 


(3.8) 


3.2 Control System Structure 

(1) Direct learning control 

In direct control, the parameters of the controller are directly adjusted to reduce the output 
error. 

(2) Indirect learning control 

system model. 
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Figure 3.2 A Neural Controller for Active Vibration Control 


Figure 3.2 depicts a basic scenario of single channel active vibration control, where F^denotes 
the transfer function (TF) of the primary channel of the system, and P 2 denotes the TT o e 
secondary channel of the system. The vibration source x , which can be measured by input 
sensors^cites the primary channel to generate uncontrollable undesirable vibration disturbance y<i 
at the output sensor location. The output vibration disturbance e is picked up by output sensorfs) 
and sent to an adaptive filter with the input sensor signal for acnve control signal generation As 
shown in the figure, the neural controller output u is fed into an actuator(secondary channe ) 2 
which produces proper anti- vibration y to cancel the vibration at output sensor location. 

For the sake of simplicity, assume that the architecture of the neural network controller is a 
single neuron with N inputs. At the k th sampling instance, the input vector consists of the N taps 
of history data of the reference signal x: [x(k), x(k-l),...,x(k-N+l)J. 

Therefore 

N-l 

u<j) = sgm[X w,xG-i)l ,3 m 

i=0 V J 

Without loss generality, let the dynamics of the secondary channel P 2 be described as: 

y(k) = P 2 [ u(k-l), u(k-2), ..., u(k-M) ] (3.10) 


OT N-l N-l -• N * ! 

y(k) = P 2 [sgm[ V wix(k-i-l)], sgm[£ w i x(k-i-2) ] sgm[ X w i x(k-i-M) ] ] 

ltS i=o 1=0 


The error signal is defined as: 
e(k) = y d (k)-y(k) 

The gradient of the performance index J = e 2 (k) is calculated as: 


(3.11) 
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(3.12) 


awk£ =2e (k)if (!0 M . 24(k) Lz<!0 

9 W| 9wi 3wi 

According to equation (3.10) 

3y(k) 3P 2 3u(k-l) t 9 P 2 9 u(k-2) | , 3 P 2 3 u(k-M) 

9 wi ” 3 u(k-l) 3 w, 3 u(k-2) 3w, 3u(k-M) 3w, 

or 

3 y(k) _y 3 P 2 3 u(k-j) 

9 wj j=i 3 u(k-j) 3 w | 

Based on the neuron model: 

= — — sgm[ X w , x(k-j-i) ] = (1- uKk-j) ) x(k-j-i) 

3 w i 9 w i i= o 

substitute this result into equation (3.13) 

8y(k) = X ( 1 - U 2(k-j) ) x(k-j-i) 

3 wj j=i 3 u(k-j) 

The update law for the proposed neural controller is then: 

W|(k+1) = Wj(k) + 2 p e(k) X ( 1 - u2 ( k 'j) ) 

j=i 3 u(k-j) 

In the case that the secondary channel can be modelled by another neuron, we have: 
M 

y(k) = sgm[ X V J u ( k ‘j)] 

J=t 

The partial derivative becomes 

— — — — = ( 1 - y 2 (k) ) vj 
9 u(k-j) 

the component of the gradient is: 

= X(I-y 2 (k))vj(l* u 2 (k-j) ) x(k-j-i) 

3w[ j=i 


i o 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


( 3 . 18 ) 


or 


d = ( 1 - y 2 (k) ) X vj ( 1 - u 2 (k-j) ) x(k-j-i) 
5wi J=i 


(3.19) 


The update law becomes: 

W|(k+l) = wi(k) 


M 


+ 2 n e(k) ( 1 - y 2 (k) ) X V J ( 1 ' u2 ( k 'J) ) x ( k ‘J _i ) 

J=i 


( 3 . 20 ) 


w. have conduct^ experimenule^uation^fthe above „°" 0 ? S 

^goriAm^o^uldplefchannei and the result will be presented in Ute next 

progress report. 
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Figure 1.1 ( f ) 





















New Actuator With Mechanical Amplifier Test 

Current Driven Mode: Voltage=0.51 V(pp), Current = 0.5 A (REM) 
Displacement = 2.2 mils (pp) 











Figure 1.2 (c) 

New Actuator With Mechanical Amplifier Test 

Current Driven Mode: Voltage=0.76 V(pp), Current = 0.75 A (REM) 
Displacement = 4 mils (pp) 















Coupling Real Time Display I Controls 

la AC /'"‘N Function Select Zoom 



New Actuator With Mechanical Amplifier Test 

Current Driven Mode: Voltage=1.0 V(pp), Current = 1.0 A (REM) 
Displacement = 5.8 mils (pp) 
















New Actuator With Mechanical Amplifier Test 

Current Driven Mode: Voltage=1.25 V(pp), Current = 1.25 A (REM) 
Displacement = 7.5 mils (pp) 














New Actuator With Mechanical Amplifier Test 

Current Driven Mode: Voltaj;e=1.49 V(pp), Current = 1.5A (REM) 
Displacement = 8.5 mils (pp) 














Current Driven Mode: Voltage=1.73 V(pp), Current = 1.75 A (REM) 
Displacement = 10 mils (pp) 
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Current Driven Mode: Voltage=1.97 V(pp), Current = 2.0 A (REM) 
Displacement =11.5 mils (pp) 



















New Actuator With Mechanical Amplifier Test 

Current Driven Mode: Vol(age=2.21 V(pp), Current = 2.25 A (REM) 
Displacement = 12 mils (pp) 













Figure 1.2 ( j ) 

New Actuator With Mechanical Amplifier Test 

Current Driven Mode: Voltage=2.45 V(pp)> Current = 2.5 A (REM) 
Displacement = 13 mils (pp) 





































